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DECAYS OF CHARMED MESONS TO PV FINAL STATES 

Bhubanjyoti Bhattacharya0 and Jonathan L. RosneJl 
Enrico Fermi Institute and Department of Physics 
University of Chicago, 5640 S. Ellis Avenue, Chicago, IL 60637 

New data on the decays of the charmed particles D°, D + , and D s to PV final 
states consisting of a light pseudoscalar meson P and a light vector meson V 
are analyzed. Following the same methods as in a previous analysis of D — > PP 
decays, one can test flavor symmetry, extract key key amplitudes, and obtain 
information on relative strong phases. Analyses are performed for Cabibbo- 
favored decays and then extended to predict properties of singly- and doubly- 
Cabibbo-suppressed processes. 

I INTRODUCTION 

In the past few years rich data on charmed particle decays have been contributed by a 
variety of experiments. Among the decays studied are those involving PV final states, 
where P and V denote light pseudoscalar and vector mesons, respectively. These decays 
obey an approximate flavor SU(3) symmetry [U [21 E], allowing one to investigate such 
questions as the strong phases of amplitudes in these decays. These strong phases can be 
important when analyzing D decay Dalitz plots in the context of studies of CP violation 
in B — > DX decays. We have recently performed a similar analysis of D —>■ PP decays [I]. 

The diagrammatic approach to flavor symmetry is reviewed briefly in Section II. Cabibbo- 
favored decays are discussed in Section III, singly- Cabibbo-suppressed decays in Section IV, 
and doubly-Cabibbo-suppressed decays in Section V. It is possible to obtain a few of the 
relevant amplitudes using factorization techniques. We discuss factorization calculations in 
Section VI and conclude in Section VII. 



II DIAGRAMMATIC AMPLITUDE EXPANSION 

A flavor-topology description of D — > PV decays uses amplitudes defined as in Ref. 
[3J. Cabibbo-favored (CF) amplitudes, proportional to the product V U( iV* s of Cabibbo- 
Kobayashi-Maskawa (CKM) factors, will be denoted by unprimed quantities; singly- Cabibbo- 
suppressed amplitudes proportional to V US V* S or V u dV* d will be denoted by primed quan- 
tities; and doubly-Cabibbo-suppressed quantities proportional to V us V* d will be denoted 
by amplitudes with a tilde. These amplitudes are in the ratio 1 : A : —A : —A 2 , where 
A = tan 0c = 0.2317 [5], with 9 C the Cabibbo angle. 
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The relevant amplitudes are labeled as T ("tree"), C ("color-suppressed"), E ("ex- 
change"), and ("A") (annihilation). For PV final states, a subscript on the amplitude 
denotes the meson (P or V) containing the spectator quark. 

The partial width T(H — > PV) for the decay of a heavy meson H may be expressed in 
terms of an invariant amplitude A as 

r(i/ - pv) = si'- 412 • (1 > 

where p* is the center-of-mass (cm.) 3-momentum of each final particle, and M# is the 
mass of the decaying particle. 



Ill CABIBBO-FAVORED DECAYS 

In Table 0] we summarize predicted and observed amplitudes for Cabibbo-favored decays of 
charmed mesons to PV. The experimental values are based on those in Ref. [5] unless noted 
otherwise. Topological amplitudes are then obtained from these processes by algebraic 
solution. The values of \Ty\ and \Ep\ are uniquely given by the rates for D s — > 7r + and 
D° —>K(j), respectively. A two-fold ambiguity then is found for the amplitude \Cp\ and 
phases of Cp and Ep, as summarized in Table HT1 

As explained in Ref. [6], the solution "B" with \Cp\ < \Ty\ is expected for a color 
suppressed amplitude. However, on the basis of fits to data from singly-Cabibbo-suppressed 
D — > PV decays, it will turn out that we will prefer the solution "A" with \Cp\ > \Ty\. In 
Fig. CD we plot these two solutions for amplitudes and relative phases of Ty, Cp and Ep. 
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Table I: Branching ratios and invariant amplitudes for Cabibbo-favored decays of charmed 
mesons to one pseudoscalar and one vector meson. 



\/l QCAT1 


i^ecay 




ID [OJ 


P 


1 A\ 




mode 




(%) 


(MeV) 


(lO" 6 ) 


D° 


K*-7T + 


T v + Ep 


5.91 ± 0.39 


710.9 


4.80 ±0.16 






Tp + E v 


10.8 ±0.7 


675.4 


7.01 ±0.23 






-jniCp — E P ) 


2.82 ±0.35 


709.3 


3.33 ±0.21 




K°p° 


"L- ( Cl\ T T^T / ^ 

V2 V V V > 


1.54 ± 0.12 


673.7 


2.66 ± 0.14 




17*° 


— — — ( C~^i T~i — 1— ~F] n /^/'t r \ 

7M\yp ' j - / v j 


96 ± 3 


579.9 


2 63 ± 41 




K T] 


-j= 6 (C P + Ep + 2E v ) 


< 0.11 


101.9 






K u> 


-±(C V + Ey) 


2.26 ± 0.4 


670.0 


3.25 ± 0.29 




— n 

TC(f> 


— Ep 


0.868 ±0.06 


520.6 


2.94 ±0.10 


D + 


m 

K 7T 


Ty + Cp 


1.83 ±0.14 


711.8 


1.68 ±0.06 




n 

TC P + 


Tp + Cy 


9.2 ±2.0 


677.0 


4.06 ±0.44 


Df 




Cp + Ay 


3.9 ±0.6 


682.4 


3.97 ±0.31 




K°K*+ 


Cy + A P 


5.3 ± 1.2 


683.2 


4.61 ±0.52 




p+7T° 


j= 2 (Ar~ A v) 




825.2 






P + V 


j=(T P -A P -A v ) 
i(2T P + A P + A v ) 


13.0 ±2.2 


723.8 


6.63 ±0.56 




p + v' 


12.2 ±2.0 


464.8 


12.5 ± 1.0 




7T + p° 


72( A v-A P ) 




824.7 






lt + L) 


-h(Ay + A P ) 


0.25 ±0.09 


821.8 


0.76 ±0.14 




7T + <f) 


Ty 


4.38 ±0.35 


711.7 


3.95 ±0.16 



Using the solutions for Ty, Cp and Ep as inputs, the other amplitudes Tp, Cy and 
Ey were obtained. The amplitude Tp was assumed real relative to Ty, in accord with the 
expectation from factorization. Six sets of solutions were obtained for each of the two cases 
| TV | < \C P \ ("A") and \T V \ > \C P \ ("B"). These solutions are listed in Table [TTT1 The 
solutions Al and A2 are found to give the best fit to the data available for singly- Cabibbo- 
suppressed D — > PV decays, and so will be singled out for special consideration. Note the 
identical magnitudes and phases of Tp, Cy and Ey in Solutions Al and Bl. 

The magnitudes and phases of solutions Al and A2 are illustrated in Fig. [2j The 
amplitudes T P + Ey = A(D° -> K~p+), Cy - Ey = V2A(D° -> K°p°), and T P + 
Cy = A(D + — > K° p + ) form a triangle whose shape is specified by their magnitudes. The 
amplitudes Cy and Ey form the sides of a quadrangle whose diagonals are Cy — Ey = 
V2A(D° -> TCp) and C v + E v = ~V2A(D° -> ICuj), and whose vertices lie on a circle 
with midpoint M. Two vertices are fixed, while the other two (A and B in Fig. [2]) lie 
at any two opposite points on the circle. An additional constraint is the magnitude of 
Cp + Ep — Ey = y/3A(D° — > K* rj). A discrete ambiguity remains, corresponding to the 
solutions listed in Tables [III and IIHI 
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Table II: Solutions in Cabibbo-favored charmed meson decays to PV final states. 







Solution A 




Solution B 


PV 


Magnitude 


Relative 


Magnitude 


Relative 


amplitude 


(io- 6 ) 


strong phase 


(io- 6 ) 


strong phase 


T v 


3.95 ±0.07 




3.95 ±0.07 




c P 

Ep 


4.88 ±0.15 
2.94 ±0.09 


S Cp t v = (-162 ±1)° 
Se p t v = (-93 ± 3)° 


2.84 ±0.09 
2.94 ±0.10 


S Cp t v = (-158.21?;°)° 
5 EpTv = (92.811;?)° 




Figure 2: Amplitudes T P , C v , and E v in solutions Al (top) and A2 (bottom). 
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Table III: Alternative solutions for Tp, Cy, and Ey amplitudes in Cabibbo-favored charmed 
meson decays to PV final states. Solutions Af - A6 correspond to \T V \ < \Cp\, while the 
solutions Bf - B6 correspond to \T V \ > \C P \ 



No. 


PV 


Magnitude 


Relative 


B(D U -> K*°r]') 




ampl. 


(f(T 6 ) 


phase 


(lO' 4 ) 


Al a 


T P 


7.46±0.2f 


Assumed 






C v 


3.46±0.f8 


5 CvTv = (172 ±3)° 


1.52 ±0.22 




Ey 


2.37±0.f9 


5 EvTv = (-110±A)° 




A2 b 


T P 


6.5f±0.23 


Assumed 






Cy 


2.47±0.22 


5 CvTp = (-174 ±4)° 


1.96 ±0.23 




Ey 


3.39±0.f6 


Se v t p = (-96 ± 3)° 




A3 


T P 


-5.67±0.22 


Assumed 






Cy 


3.64±0.27 


Sc v t p = (-46 ± 4)° 


1.42 ±0.28 




Ey 


2.09±0.28 


Se v t p = (-122l|)° 




A4 


T P 


-5.60±0.24 


Assumed 






Cy 


f .68±0.24 


&c v T P = (-20 ± 6)° 


2.21 ± 0.25 




Ey 


3.85±0.f5 


Se v t p = (-94 ± 3)° 




A5 


T P 


-3.22±0.2f 


Assumed 






Cy 


f .79±0.32 


Sc v t p = (-104 ±5)° 


2.18 ± 0.25 




Ey 


3.79±0.f3 


Se v t p = (-180l|)° 




A6 


T P 


3.2f±0.2f 


Assumed 






Cy 


f .78±0.3f 


5 Cv t p = (105 ± 5)° 


2.18 ± 0.25 




Ey 


3.80±0.f3 


Se v t p = (-18011)° 




Bf 


T P 


7.46±0.2f 


Assumed 






Cy 


3.46±0.f7 


5 CvTp = (172±3)° 


0.33±0.05 




Ey 


2.37±0.f9 


5 Ev t p = (-110 ±4)° 




B2 


T P 


6.43±0.22 


Assumed 






Cy 


3.95±0.24 


5 CvTp = (-143±4)° 


u.uoz_ 021 




Ey 


f.40±0.32 


5 EvTp = (-711?)° 




B3 


T P 


4.53±0.24 


Assumed 






Cy 


0.80±0.2f 


r ( 1 o n + 1 6 \ o 

5 Cv t p = (1301^)° 


I.l8±0.l0 




Ey 


4.12±0.15 


5 EvTp = (72 ± 3)° 




B4 


T P 


4.97±0.22 


Assumed 






Cy 


3.28±0.29 


5 CvTp = (126 ± 4)° 


0.42±0.l0 




Ey 


2.6f±0.25 


6 EvTp = (47±5)° 




B5 


T P 


-3.33±0.22 


Assumed 






Cy 


0.75±0.f9 


5 CvTv = (16411*)° 


l.l9±0.ll 




Ey 


4.13±0.17 


S EvTv = (-140 ± 2)° 




B6 


T P 


-7.70±0.21 


Assumed 






Cy 

Ey 


4.01±0.17 
f.24±0.22 


S Cv t v = (17l|)° 
8e v t v = (-52lg)° 


0.020±0.0ll 



"Preferred solution based on fit to singly-Cabibbo-suppressed decays. 

6 Alternative solution giving acceptable fit to singly-Cabibbo-suppressed decays. 
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Predictions for the branching ratio for D° — ► K rj', listed in the last column of Table 
IIVt in principle allow one to distinguish among various solutions. In addition, we shall 
see that only solutions Al and A2 give rise to acceptable fits to singly-Cabibbo-suppressed 
decays. 

We now state a relationship between \Tp\ and Cabibbo-favored D s decay amplitudes: 

\A(D S - p + r/)| 2 = \T P \ 2 + \A{D S -> n + cu)\ 2 - \A(D S - p + V )\ 2 (2) 

Using the value of \Tp\ from solution Al of Table II I II and the decay amplitudes (D s — ► 
p + ?7, 7r + u;) from Tablefl], we calculate the amplitude: \A(D S — > p + r)')\ = (3.50±1.15) x 10~ 6 , 
which deviates from the experimental value (Table [T|) by a large amount. This problem 
with the quoted experimental rate for D s — > p + rj' was already noted in Ref. [6]. It indicates 
either the importance of neglected amplitudes involving the flavor-singlet component of rf , 
or an overestimate of the experimental decay rate in this mode. 

The remaining parameters Ap and Ay were determined using the amplitudes of D s — > 
(iT K+, K K* + , it + u) and have been listed in Table IIVI A direct calculation of the 
amplitudes for D s — > p + (r], r]') is now possible using these amplitudes. For the ampli- 
tude solutions (Al, A2) preferred by fits to singly-Cabibbo-suppressed decays, we find 
B(D S — > p + rj) = (5.6 ± 1.2, 5.55 ± 0.60)%, to be compared with the experimental value of 
(6.63 ± 0.56)%, and B(D S -> p+rf) = (2.9 ± 0.3, 1.89 ± 0.20)%, to be compared with the 
experimental value of (12.5 ± 1.0)%. The agreement between prediction and experiment 
for B(D S — > p + i]) is good for the solutions Al, A2, Bl, and B2, while no solution gives 
agreement for B(D S — > p + f]'). We await forthcoming CLEO data on this mode. 

IV SINGLY-CABIBBO-SUPPRESSED DECAYS 

The topological amplitude decomposition of singly-Cabibbo-suppressed decays of D° — > PV 
is listed in Table |V] along with the measured branching ratios and amplitudes for the 
decays. Unlike the D — > PP case [1], here we have neglected the Okubo-Zweig-Iizuka 
(OZI) suppressed disconnected diagrams that form the Singlet-Exchange {SE') and Singlet- 
Annihilation (SA') amplitudes. 

We now make use of the amplitudes determined in Section III to predict the singly- 
Cabibbo-suppressed decay amplitudes. Here we assume the simple hierarchy of amplitudes 
explained in Section II. Based on the available data we calculated the global x 2 °f singly- 
Cabibbo-suppressed D — > PV decays for solutions A1-A6 and B1-B6. Solutions Al and 
A2 have the two lowest values of x 2 an d hence were chosen as the preferred and alternative 
solutions. Table IVT1 summarizes the global x 2 values for each of the twelve solutions. It also 
includes, for each solution, two processes that contribute the most towards a high value of 

x 2 . 

One notes in Table I VI I that the main processes contributing to high global x 2 f° r an 
solutions are D° — > <J)tt and D° — > p°7r°. The solutions Bl-6, which correspond to \Cp\ < 
\Ty\, yield high x 2 f° r th e process D° — > 07r°. The amplitude of this process depends 
only on Cp' . This shows that \Cp\ < \Ty\ is not favored by the process D° — ► 07r°. The 
processes D° — ► p°7r° and D + — > p°7i + contribute to high x 2 f° r the solutions A3-6. 

The predicted and experimental D° branching ratios are in qualitative agreement but 
with some notable exceptions. The predictions for D° — > irp fall slightly short of experiment 
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Table IV: Solution for annihilation amplitudes in Cabibbo-favored charmed meson decays 
to PV final states. 



No. 



PV 
amplitude 



Magnitude 

(io- 6 ) 



Relative 
phase 



T36™ 

1 9 r +0.34 
i - zo -0.31 



-i5ii^)° 



Prediction 

(%) 



Al a 



A P 
A v 



5a p 



A9ll°) 



-0.22 
-0.18 
i 9 o+0.31 



B(Df 



5.6 ± 1.2 
2.9 ±0.3 



A P 
A v 



2.15; 



5a p 

^4v 



-1791 9 , 
-19l?|)° 



o.96±8;g 



-14j 
+21\o 

u) 



B(D; 
B(Dj 



rj' p+ 



5.55 ±0.60 
1.89 ±0.20 



A3 



A P 
A v 



$A X 



= -89 



-0.21 
q ori+0.23 
°- zu -0.19 



(34 



-109±f)° 



B(Df 



Tj' p H 



4.20 ±0.81 
1.45 ±0.28 



A4 



A P 
A v 



4.271 



8a v 



-72 



+6\o 



1 QC+0.35 
^•oo_n 



(-123^)° 



B(Df 



> r]p^ 



2.77 ±0.27 
1.77 ±0.18 



A5 



A P 
A v 



1.93 



-0.24 
hl.21 
-0.27 



+0.22 

-0.31 

2 85+ - 21 

Z.OO_ 26 



5 Av = (Q9tfy 



B(Df 
B(D+ 



• T) p^ 

v'p" 



0.58 ±0.06 
0.70 ±0.07 



(-36 ±7)° 



A6 



A P 
A v 



2. 



5a 
5a, 



B(Df 



> 7]p^ 



1.61 ±0.17 
0.43 ± 0.04 



-0.32 
, 74 +0.44 
l - ' ^-0.28 



:±i2i^ y ) c 

-96^)° 



Bl 



A v 



1.57! 



$A P 



1 qpr+0.51 
l.OO_ 27 

1 52+ - 70 

i - dz -0.21 



(-74lf )° 
(±1501^)° 



B(£>+ 



> rip-" 

Tj' p^ 



7.08 ± 1.03 
2.53 ±0.37 



O.Oo_2 ii 

1 Sfi+ - 70 
J-SO-O^ 



B2 



Ay 



0.24 
37 
22 



B(Df 
B(D+ - 



' VP 



±mif)° 



B3 



A P 
A v 



3.85 
2.7818 



( 



i 74 +o.34 

x - ' 4 -0.23 

1.16j 



= (+77^5)° " 
;-140 ± 12)° 



B(D1 
B(Df 



> r)p^ 
i P^ 



2.42 ±0.16 
1.01 ±0.07 



B4 



A P 
A v 



$A f 



-0.27 
-0.23 



^4, 



o OO+0.29 
°- zz -0.38 



(-eolfO 



B(Df - 



77' p" 1 



3.04 ±0.70 
1.18 ±0.27 



B5 



Ap 
A v 



<$Af 
$A X 



B(£>+ - 



► ?7P 
77' p^ 



+ 



u -° ' -0.29 
1 90+O.23 



: (+45^)° 
(+1681^)° 



1.30 ±0.10 

o.57ii8:8g 



B6 



A P 
A v 



5a 
5a 



p 

V 



B(Df 
B(Df - 



v' P^ 



4.80 ±2.54 
3.42 ± 1.81 



"Preferred solution based on fit to singly-Cabibbo-suppressed decays. 

b Alternative solution giving acceptable fit to singly-Cabibbo-suppressed decays. 
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Table V: Branching ratios and invariant amplitudes for singly-Cabibbo-suppressed decays 
of charmed mesons to one pseudoscalar and one vector meson. 



Meson 


Decay 


Representation 


B [5J 


V 


1 A 1 




moQe 






l^ivie v j 


l iU ) 


no 




(rp 1 , 771 i\ 


U.4y / rrU.U/o 


I Oo.o 


i.ZDrrU.Uo 




A* 


_(7V + Ev'\ 

\ J- p ~ J-zy ) 


Q80±0 040 


763.8 


1 76±0 04 




7T /) 

r 


— ( h*l T 1 ' — 1 — A'/'t / ' C~^. 1 i ' ( / ' ^ 

2 P i ^ r J 


373±0 022 


764.2 


1 08±0 03 






Ty ' + Ep ' 


0.153±0.015 


609.8 


0.97±0.05 




K* + 


Tp' + Ey' 


0.441±0.021 


609.8 


1.65±0.04 




K TT° 


Ey ' — Ep 


< 0.18 


605.3 






K°K*° 


Ep ' — Ey 


< 0.09 


605.3 








l ri i 
~7n L, p 

iii i i 

2 {Ep ' + Ey ' — Cp' + Cy ') 


0.124±0.012 


644.7 


0.81±0.04 








761.2 






7] p° 


-^{2Cy' - C P ' - Ep' - Ey') 




652.0 






r\ijj 


-±{2Cy' + C P ' + E P ' + Ey') 




488.8 






7](j) 


-75 (Cp' — Ep' — Ey') 




648.1 






r)'p° 


^{Ep' + Ey' + Cp' + Cy') 




342.5 






rj'uj 


{Ep ' + Ey ' + Cp' — Cy') 




333.5 




D + 


p°7T + 


j-(A P ' -Ay' -C P ' -Ty') 


0.082±0.015 


767 


0.32±0.03 




UJ 7l + 


-l(A P ' + Ay' + Cp' + Ty') 


< 0.034 


764 






(f) 71" + 


Cp' 


0.620±0.070 


647 


1.13±0.06 




K*°K + 


{Ty' -Ay') 


0.435±0.048 


611 


1.03±0.06 






MAy' -A P ' -Cy' -T P ') 

^{A v ' + Ap' + 2Cy' + Tp') 




767 








< 0.7 


656 






v'p + 


j= 6 {Cy'-Ay'-Ap'-T P ') 


< 0.5 


349 






K°K* + 


{Tp'-Ap') 


3.18±1.38 


612 


2.78±0.60 


Dt 


7T° K* + 


{Ay'-Ty') 
--L(Cy' + Ay') 


0.225±0.039 


773 
775 


0.79±0.07 




7]K* + 


MT P ' + 2Cy' + A P ' -Ay') 
-±{2T P ' + Cy' + 2A P ' + Ay') 




661 






rj'K*+ 




337 






K°p+ 


{Ap'-Tp') 




743 






K+p° 


-j= 2 {C P ' + Ap') 


0.27±0.05 


745 


0.92±0.09 




K+u 


-A_(C P '-A P ') 

Ty' + Cp' + Ay' 




741 






K + (j) 


< 0.057 


607 
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Table VI: Global x 2 values for fits to singly-Cabibbo-suppressed D — > PV decays. Also 
included are the process that contribute the most to a high x 2 value. 



No. Global Worst Processes (High Ax 2 value) 

X 2 Decay Channel B th {%) B expt {%) A X 2 

Al" 6L8 D + -»• K* u K+ 0.17 ±0.04 ~ 0.435 ± 0.048 16.1 
D + -> m 7T+ 0.16 ±0.04 < 0.034 11.6 

A2^ 65^9 D+ -»• 7T U K+ 0.17 ±0.03 0.435 ± 0.048 21.4 

D° -> p° tt° 0.27 ±0.02 0.373 ±0.022 10.1 

A3 341.4 D° -»■ p u 7T° (4.3 ±3.1) x 10~ 3 0.373 ±0.022 275.2 

£> + ^p°7r + (1.5 ±4.0) x 10~ 3 0.082 ±0.015 25.1 

A4 1671 D° -»• p u 7T° 0.12 ±0.01 0.373 ± 0.022 95~X 

£>+ -> p° 7T+ 0.73 ±0.12 0.082 ±0.015 31.4 

A5 3241 D () -»• p u tt u (6.1 ±3.1) x 10" 3 0.373 ±0.022 272.6 

D + ^ K*°K° 0.19 ±0.02 < 0.09 11.9 

A6 149.8 D+ -»■ p° 7T+ 0.91 ±0.09 0.082 ±0.015 51.1 
D + ^K*°K + 0.12 ±0.03 0.435 ±0.048 32.1 

Bl 244.0 L>° -»• p u 7T° 0.12 ±0.01 0.373 ±0.022 95.3 

£>° ->■ 07T° 0.042 ±0.003 0.124 ±0.012 45.3 

B2 1551 D° -»• 07T° 0.042 ± 0.003 0.124 ±0.012 45X 

£>+ -> 07T+ 0.21 ±0.01 0.62 ±0.07 32.9 

B3 165.7 L>° ->• 07T° 0.042 ± 0.003 0.124 ±0.012 45.3 

£>+ -» 07T+ 0.21 ±0.01 0.62 ±0.07 32.9 

B4 15L7 L>° -»• 07T° 0.042 ± 0.002 0.124 ±0.012 45X 

D+ -> p° 7T+ 1.44 ±0.23 0.082 ±0.015 34.4 

B5 518.8 £>° -»• p u 7T° (5.4 ±2.8) x 10~ 3 0.373 ±0.022 274.8 

D + -> p° 7T+ 1.71 ±0.21 0.082 ±0.015 59.3 

B6 401.3 D° -> p u 7t° 0.015 ±0.006 0.373 ± 0.022 245.9 

D°^07r 0.042 ±0.003 0.124 ±0.012 45.3 



"Preferred solution. fc Alternative solution. 
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for all charge states, most prominently for 7r°p°. Recall that the predicted branching ratio 
for D° — > 7r + 7T _ lies significantly above the experimental value |4J. The predictions for 
D° -> A+A*~ and D° -> K~ K* + are not badly obeyed, while those for D° -> A^A* and 
.D — > X A are far below the current experimental upper limits. The predicted branching 
ratio for D° — > 7T°(/> is approximately same as the observed value. The value of x 2 f° r 
solutions # Al and A2 are respectively 61.8 and 65.9 (Table fVT|) . where we have used the 
18 data points for which the branching ratios are available. 

In Table IVHl we present our predictions for branching ratios of singly- Cabibbo-suppressed 
D — > PV modes corresponding to the two solutions Al and A2 having the lowest value of 
global x 2 f° r these modes. There is little one can do to distinguish between them given 
the available data on branching ratios. Both solutions yield fairly similar central values for 
most of the singly- Cabibbo-suppressed D —>■ PV modes. A slight distinction may be made 
in a few cases. For example, the predicted central values of B(D° — > (A A*°, A° A* )) are 
larger for solution Al than for A2, though differing only by 1.5a. Another example is the 
process D° — > tt°u, for which the central value of the branching ratio in solution A2 is nearly 
three times its value in Al. Still another example is the process D + — > i]'p + , for which 
the predicted (very small) branching ratio in Al is twice its value in A2. Measurements of 
the branching ratios for both Cabibbo-favored and singly-Cabibbo-suppressed decays with 
higher precision will be necessary in order to distinguish between the two solutions. 

V DOUBLY-CABIBBO-SUPPRESSED DECAYS 

We now characterize the doubly-Cabibbo-suppressed or wrong-sign (WS) decays of D — ► 
PV. A detailed list of possible decays and the corresponding topological amplitude de- 
compositions are given in Table IVIII1 We used the Cabibbo-favored amplitudes calculated 
in section III to predict the WS amplitudes, using the simple hierarchy of amplitudes as 
explained in Section II. The predicted amplitudes have been included in Table rvTTTI for the 
preferred (Al) and alternative (A2) solutions. 

The experimental values for the following decays are available in the literature [5]: 

B(D° ^ K*+7T~) = (3.0i?;!j) x 1(T 4 (3) 
B{D + -f K*°7i+) = (4.35 ± 0.9) x 10~ 4 (4) 

The predicted values for these branching ratios (Table IVIIII) are in satisfactory agreement 
with the experimental values quoted above. An interesting point to note is that both 
solutions Al and A2 give the same predicted central values for these branching ratios, but 
A2 has a larger error bar on both of them. Several other branching ratios in Table IVIIII 
predicted to exceed 10~ 4 may help to distinguish between solutions Al and A2. These 
include B(D° -> A*°7r°), B(D+ -> K* + n°), and B(D + -> K + p°). Reduction in errors 
on predictions will be needed in order that these distinctions exceed 2-2. 5a. Some of the 
doubly-Cabibbo-suppressed decays in Table IVIIII may be observable in Dalitz plots of D 
decays to three pseudoscalars through interference with Cabibbo-favored PV decays. For 
example, D° — > As7r + 7r~ might be able to provide new information about the decay process 
D° — > K* + 7i~ , while D + As7r + 7r° could provide information about D + — > K* + tt°. 
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Table VII: Comparison between predicted amplitudes based on Cabibbo-favored decays and 
the experimental values for singly-Cabibbo-suppressed decays of D° to a pseudoscalar and 
a vector meson. Predictions are listed for preferred (Al) and alternative (A2) solutions. 





PV Decay 


Experimental 


Predicted B (%) 


Meson 


Mode 


£ (%) 


Solution A1 


Solution A2 


D° 


H 


D 497 ± f) 023 


U.oy ± U.Uo 


n qo -i- n nQ 
U.oy it U.Uo 




H 




u.o^ it u.uu 


n ra 4- n nfi 

u.o^ it u.uu 




7T° p° 
r 


0.373 ± 0.022 


29 ± 02 


27 ± 02 






0.153 ±0.015 


20 ± 01 


20 ± 01 






0.441 ±0.021 


43 ± 03 

v ./ . 7t (J _1_ V . V ".J 


43 + 03 






< 0.18 


u.uuou zn u.uuou 


n nn9n 4- n nm <s 

u.uuou in u.uuiu 




IV IV 


<r n HQ 
- - u.uy 


u.uuou zn u.uuou 


n nn9n 4- n nm r 

u.uuou in u.uuiu 




7T°0 


0.124 ±0.012 


U.LZZ It U.UU I 


n i oo -i- n nn7 

U. 1ZZ It U.UU ( 




7T°CJ 




U.U4o it U.UUo 


n 1 1 q 4- n ni o 
u.iiy it u.uiz 




'/ I 




U. 1UO It U.Ulo 


u.uyo it u.uiu 




Tj UJ 




n i /in 4- n nriQ 
u. i4u it u.uuy 


n 1 97 4- n nna 
u. lz i it u.uuy 




71 <b 
'/ Y 




n no.'? -i- n nno. 
u.uyo it u.uuy 


n i a 4- n m 

U. 14 It U.Ul 




'/ r 




n m ^zl 4- n nnno, 
u.uio4 it u.uuuy 


n ni 4- n nnno 
u.uioo it u.uuuy 




Tj ' UJ 




U.UUOO it u.uuuo 


n nn77 4- n nnnc; 
U.UUrf ( it u.uuuo 


D + 




082 ± 015 


0.097 ±0.048 


0.23 ±0.12 




C<J 7T + 


< 0.034 


0.15 ±0.04 


0.14 ±0.12 




71" + 


0.620 ±0.070 


0.62 ±0.04 


0.62 ±0.04 




K*°K+ 


0.435 ± 0.048 


0.17 ±0.04 


0.17 ±0.03 




ttV 




0.062 ±0.047 


0.012 ±0.015 






< 0.7 


0.0017 ±0.0040 


0.0057 ±0.013 




Vp+ 


< 0.5 


0.083 ±0.010 


0.044 ±0.005 




K°K* + 


3.18 ±1.38 


1.66 ±0.20 


1.66 ±0.12 


Df 




0.225 ±0.039 


0.15 ±0.04 


0.15 ±0.03 




TT°K* + 




0.049 ±0.012 


0.020 ±0.008 




T]K*+ 




0.014 ±0.011 


0.012 ±0.008 




r]'K* + 




0.029 ±0.006 


0.015 ±0.003 




K°p+ 




1.29 ±0.15 


1.29 ±0.09 




K+p° 


0.27 ±0.05 


0.33 ±0.05 


0.42 ±0.05 




K + uo 




0.108 ±0.029 


0.072 ±0.033 






< 0.057 


0.038 ±0.009 


0.037 ±0.028 
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Table VIII: Branching ratios and invariant amplitudes for doubly- Cabibbo-suppressed de- 
cays of charmed mesons to one pseudoscalar and one vector meson. Predictions are shown 



for favored (Al) and alternative (A2) solutions. 


Meson 


Decay 


Representation 


P* 


Predicted B{1Q- A ) 




mode 




(MeV) 


Solution Al 


Solution A2 


D° 


K*+ 7T- 


Tp" + E v " 


711 


3.63 ±0.26 


3.63 ±0.27 




K*° 7T° 


(C P "-E V ")/V2 


709 


0.55 ±0.06 


0.80 ±0.08 




K*°r] 


(C P " -E P " + E V ")/V3 


580 


0.38 ±0.04 


0.37 ±0.04 




K*°r)' 


-(C P " + 2Ep" + E v ")/V6 


102 


0.0046 ±0.0004 


0.0052 ±0.0004 




K+ p- 


T v " + E P " 


675 


1.46 ±0.10 


1.46 ±0.10 




K°p° 


(C V "-E P ")/V2 


674 


0.70 ±0.07 


0.39 ±0.05 




K°u 


-(C v " + Ep")/V2 


670 


0.58 ±0.06 


0.52 ±0.06 




K°<f> 


—Ey " 


521 


0.16 ±0.03 


0.33 ±0.03 


D+ 


K*°7I + 


C P " + Ay" 


712 


2.94 ±0.52 


2.94 ±0.65 




K*+ 7T° 


(T P "-A V ")/V2 


714 


3.74 ±0.49 


2.71 ±0.30 




K* + r] 


-(Tp"-A P " + Ay")/VS 


586 


3.37 ±0.43 


3.37 ±0.25 




K* + T]' 


(Tp" + 2A P " + Ay")/V6 


137 


0.0095 ±0.0029 


0.0026 ±0.0010 




K°p + 


Cy" + A P " 


677 


3.43 ±0.75 


3.43 ± 0.47 




K+p° 


(Ty"-A P ")/V2 


679 


2.17±0.40 


3.01 ±0.24 




K+uj 


(Ty" + A P ")/V2 


675 


0.64 ±0.21 


0.26 ±0.07 




K + <P 


Ay" 


527 


0.12 ±0.06 


0.12 ±0.06 


Df 


K*°K + 


Ty"+C P " 


682 


0.20 ±0.02 


0.20 ±0.02 


K*+ K° 


T P "+Cy" 


683 


1.18 ±0.16 


1.18 ±0.18 
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VI FACTORIZATION COMPARISONS 



In the current section we compare our results for the amplitudes of Tp and Ty with the 
values extracted from explicit evaluation of the tree diagram assuming factorization [Jj. In 
order to calculate Tp we use the decay D° — > K~ p + . In this scenario the spectator u quark 
goes from D° to the pseudoscalar K~ and so we use the standard form of the (D — > P) 
current [8]: 

H P = f+(q 2 )(pD +Pk)^ ~ f-(q 2 )(pD-pK)^ (5) 
where / + and /_ are the relevant form factors. The current we use for the vector meson is 

= e»m p f p (6) 

where e M represents the polarization of the vector meson, m p is its mass and f p is the 
associated decay constant. The invariant amplitude and the decay rate for the process 
D° — > K~ p + via the tree diagram may then be written as 

M(D°^K-p + ) Tp = (7) 

T{D^K-p+) Tp = £ \M{D«^ K~ p+) Tp \ 2 (8) 

After summing over the p polarization and taking the modulus squared of the invariant 
amplitude one obtains the final form for \Tp\: 

lrr , G F \V ud \\V cs \\U(m 2 p )\f P 



y/2 P* 



(9) 

(m 2 D — m? K ) 2 — m 2 p (m 2 D + m? K + 2m D \Jrn? K + p* 2 ) (10) 

= (5.45 ± 0.07) x 1(T 6 (11) 

which is to be compared with the values quoted in Table IHIl and favors solution A2 over 
Al. 

In obtaining the result stated above we used |/ + (mp)||V^a| = 0.869 ± 0.009 [12]. The 
particle masses and the quantity jV^I were taken from [5]. p* is as quoted in Table fl] We 
calculated the value of f p using the following formula: 



fp fir 



B(t — > v T p 



in- - in- ni T 



m 2 — m 2 



m 2 + 2m 2 



= (209±1.6)MeV (13) 

where once again the particle masses and branching fractions were taken from [5]. 

A similar approach may be taken in order to evaluate \Ty\ by looking at the decay 
D° — > K*~ 7r + via the tree diagram. In this case the spectator u quark goes from D° to the 
vector meson K*~, so we use the standard forms of the (D — > V) vector and axial- vector 
currents [8] and the pion current [9]: 

V p = ige ppTa e' p (p D + p K ,y(p D -p K *y (14) 
A p = f e p, + a +( e *-PD)(PD+PK*)p + a~(e*-PD)(PD-PK*)t, (15) 
^ = (16) 
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We obtain for the amplitude \Ty\ the following expression: 

\ T v\ = ^\Vcs\\V ud \U—^-\f + a + (m 2 D - m 2 K ,) + a_m^| (17) 
V2 rn K * 

In principle this can be used to calculate Ty once the form factors are given. However, we 
may adopt a simplification using a result from Ref. [6], based on the earlier discussion in 
Ref. [9]. In the heavy-quark limit one expects T(D — > A'*7r + )j> = r(P — > Ktt + )t and hence 

^|T|L = (p* K ^f\T\ 2 K ^ , (18) 

where 7V*7r = Ty. In Ref. [1] we found in a fit to D — > PP amplitudes that (Tl^ = 
(2.78 ± 0.13) x 1(T 6 GeV. With p* Kw = 0.861 GeV and p* K ^ = 0.711 GeV we then obtain 
the result 

T v = (4.3 ±0.2) x 10" 6 , (19) 

in reasonable agreement with the value of (3.95 ±0.07) x 10~ 6 quoted in Table HTl especially 
considering the uncertainties associated with QCD corrections and with the use of the 
heavy-quark limit for the final strange quark. 



VII CONCLUSIONS 

We have used the flavor topology description to study the validity of flavor SU(3) for de- 
scribing D — > PV decays, to obtain relative phases and magnitudes of various contributing 
amplitudes, and to predict rates for as-yet-unseen singly- and doubly- Cabibbo-suppressed 
decays. We assumed flavor £77(3) to be an exact symmetry for the tree level diagrams. 
We found that singly- Cabibbo-suppressed decays favor a ratio of color-suppressed to tree 
amplitudes \Cp/Ty\ > 1, where the subscript denotes the meson (P or V) containing the 
spectator quark. The present data for the Cabibbo-favored decays are compatible with 
twelve distinct sets of solutions for the amplitudes T P , Cp, E P , A P , T v , Cy, E v , and 
Ay (up to discrete ambiguities). However, on the basis of experimental branching ratios 
for singly-Cabibbo-suppressed decays we were able to choose two sets of solutions giving 
substantially lower values for \ 2 than the other ten. 

Our predictions of the branching ratios for singly-Cabibbo-suppressed decays deviate 
from the available experimental data in several cases, such as those in the first four lines of 
Table VI. This shows that flavor £77(3) is not an exact symmetry. However flavor ££7(3) 
breaking, though present, is no worse in D — > PV decays than in the D —>■ PP decays 
discussed in Ref. jl]. 

Our prediction for the —>■ rj' p + branching ratio is much lower than the available 
experimental value. Either there are additional contributions to rf production which we 
have neglected, or the experimental situation needs to be re-evaluated. 

Our analysis of the singly- Cabibbo suppressed decays shows that processes such as 
D° — > n°u can be used to distinguish between the two most likely amplitude solutions. 
The mean values predicted for the branching ratios of these processes differ by nearly a 
factor of three in the two solutions, but experimental data are not yet available to resolve 
this problem. 
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The branching ratios predicted for doubly-Cabibbo-suppressed decays are close to the 
experimental values in the two cases for which data are available. A precise measurement of 
a few of the other branching ratios may help select one of the two most-favored amplitude 
solutions. 

Finally, factorization computations of the tree amplitudes agree with results obtained 
in direct analyses. However, a more precise calculation of the amplitudes using the factor- 
ization assumption could be done if data on the relevant form factors were available. 
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